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O ! Abstract 

cn . 

Starting from the general Lindhard theory describing the partition of particles en- 
ergy in materials between ionisation and displacements, analytical approximate so- 
^ ■ lutions have been derived, for media containing one and more atomic species, for 

C^- , particles identical and different to the medium ones. Particular cases, and the limits 

rr \ of these equations at very high energies are discussed. 

O. 
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1 Introduction 



The characteristics of semiconductor devices and other crystalline materials 
used in high fluences of particles are strongly affected by the effects of radia- 
tion. In the recent years, important results have been achieved in the radiation 
damage analysis, but not all the degradation mechanisms have been com- 
pletely understood up to now. After the interaction of the radiation field with 
the semiconductor material, mainly two classes of degradation effects were 
observed: surface damage and bulk material damage, due to the displacement 
of atoms from their sites in the lattice. 
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If heavy particles as pions, protons, neutrons, ions, etc., produce both sur- 
face and bulk damages, for electrons and gammas the effects are dominantly 
produced at the surface. The surface effects are not a main obstacle to semicon- 
ductor device operation, while the bulk ones have a much greater importance. 
The disruption of the symmetry of the crystal, consequence of the bulk effects, 
causes the formation of energy levels in the normally forbidden region between 
the valence and conduction bands, altering the material properties. 

The process of partitioning the energy of the recoil nuclei (produced due the 
interaction of the incident particle with the lattice site nucleus) in new inter- 
actions processes, between electrons (ionisation) and atomic motion (displace- 
ments) is considered in the general form, in the frame of the Lindhard theory 

[!]■ 

Up to now, there exist in the literature some calculations of the Lindhard 
curves for some crystalline materials: diamond [2,3], Si [4,5], GaAs [6], InP 
[7], and GaP, InAs, InSb [8], in Al and some scintillator materials [9]. 

The knowledge of the Lindhard energy partitioning curves is essential, together 
with the detailed nuclear interaction mechanisms, in the correct evaluation of 
the concentration of primary defects produced in materials that work in intense 
fields of radiation. In this paper, we present the concrete equations necessary 
in the analytical calculations of these curves. 



General hypothesis and analytical approximations of Lindhard 
equations 



The general process considered in the study of the interaction between the 
incident particle and the solid is the following: the particle, heavier than the 
electron, with electrical charge or not, interacts with the electrons and with 
the nuclei of the semiconductor lattice. It loses its energy in several processes, 
which depend on the nature of the particle and on its energy. The effect 
of the interaction of the incident particle with the electrons of the target is 
the ionisation. The quantity characteristic for this process is the energy loss 
dE/dx (stopping power). The nuclear interaction between the incident particle 
and the lattice nuclei produces bulk defects. As a result of the interaction, 
depending on the energy and on the nature of the incident particle, one or more 
light particles are formed, and usually one (or more) heavy recoil nuclei. The 
nucleus has charge and mass numbers lower or equal with that of the medium. 
After this interaction process, the recoil nucleus or nuclei are displaced from 
the lattice positions in interstitials. Then, the primary knock-on nucleus, if 
its energy is large enough, can produce the displacement of a new nucleus, 
and the process continues as long as the energy of the colliding nucleus is 



higher than the threshold for atomic displacements. This phenomenon can 
be regarded as a cascade process. We denote by primary displacements all 
the displacements produced as a result of primary interactions, without any 
further rearrangement of the vacancies and interstitials. The physical quantity 
characterising the process is the concentration of primary defects (or related 
quantities, for example the non ionising energy loss) produced per unit of 
fluence of the incident particles. 

In all the subsequent discussions, the primary recoil will be the particle who's 
energy partition is to be calculated. As specified before, the primary recoil is 
either a nucleus of the medium, or a nucleus with a lower mass and charge 
numbers. As a consequence, for each medium a whole family of curves can be 
obtained. Also these Lindhard curves can be directly used in the evaluation 
of the damage produced in materials by ion beams, if the energy of particles 
in the beam is identified with recoils energy. 

The incident particle has an initial energy E, and, due to the interaction with 
the target during the slowing down, this energy is transferred, on the one side 
to atoms (the quantity Ei), and on the other side to electrons (the quantity 
E2). Obviously: 

E = Ei + E2 (1) 



In accord with [1], the equation satisfied by Ei, as a function of the energy E, 
is: 






(2) 



where dan^e are the differential cross sections corresponding to particle scat- 
tering on nucleus and electrons, T„ (Tgj) represent the energies transferred to 
nucleus (and respectively to electrons), U is the energy wasted in disrupting 
the atomic binding, and f/j are the corresponding ionisation energies, with 
summation over all electrons. 

The equation takes into consideration all the steps of the interaction until the 
transferred energy is lower than the threshold energy for displacements. In 
this general form, the equation is practically impossible to be solved. It has 
been simplified [1], using some physical approximations: the electrons do not 
produce recoil nuclei with appreciable energy, so that the function Eie (for 
electrons), can be obtained separately; the binding energy of the atom in the 
lattice can be neglected; the energy transferred to electrons and respectively 



to nuclei is small in respect to the particle energy, and electronic and nuclear 
collisions can be separated. More, the energy transferred from electrons to 
nuclei is negligible, Eie = 0. 

Using the up - mentioned approximations, the equation for the function Ei 
becomes: 



{Sn + S,)E[= j E,{T)^dT (3) 



where Sn,e = J Tn,ed<7n,e represent the stopping cross sections for nuclei and 
electrons. The boundary condition for the equation 3 is: Ei (E) / E — * 1 for 
^^0 . 



The parameter ^ (_E), defined as: 



i{E) = ^ (4) 



represents a measure of the division of energy dissipation into electronic and 
atomic motions. For the simplest case, of a medium consisting of only one 
atomic species, and of the particle - primary recoil identical to the particles of 
the medium ( Zpart = Zmed = Z, and Apart = Amed = A), and for an electronic 
cross section Se oc E'^^'^ and with S'„ derived from a Thomas-Fermi poten- 
tial, it was shown [1] that for C, (E) there exist roughly three distinct energy 
regions. In the first one, the nuclear stopping is dominating, and relatively 
little energy goes into electronic motion; in the second region, the nuclear 
stopping start decreasing, while the electronic one increases as E^^"^, so the 
quantity C, (E) increases rapidly, and the fraction of the energy that goes into 
electronic motion increases correspondingly; in the third region, the electronic 
stopping starts decreasing, while C, {E), though still increasing, approaches an 
asymptotic value. 

It is convenient to look for analytical approximate solutions. The simplest case 
to treat mathematically is that of a power low potential, V (r) oc r"'', corre- 
sponding to nuclear scattering. In this case, the equation 2 can be replaced 
with a differential equation whose solution is the hypergeometric function. In 
the first region an potential represents a good approximation. The solution is: 

E, (E) oc E, {-12 + 6 [l + 2 (EJE^^] ■ log (l + (EJE^^) } (5) 



where Ec cc Z ■ A is the upper limit of the first energy region. In the limit 
E/Ec ^ 1, the solution for Ei (E) can be expressed even simpler, as a power 
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series in (E/Ec) . The solution described by equation 5 is a fast increasing 
one on energy. 

In the second region, the cross section could be modelled using the Rutherford 
scattering formula, i.e. using s = 4/3. The analytical formula derived for Ei 
[1] is: 



il/l £/ OC Ci-C/fe < 1 ^f ' log ; -1= ■ 

^ ' \ 4v^^ ^ ei^-^i^^ + i 



^ ^y^^rct^ni^^) +C,r"^E~''^ (6) 



where ^ {E) defined by eq. 4 is given by E/Ei,, E/, being the energy at which 
the two stopping cross sections are equal; Ef, has the same dependence on A 
and Z as E^, Ci is of the order of unity, and has the approximate expression: 

12 (a; + l)(a;-2) , , 
Ci ~ — + 6^ 4 '- ■ log {x + 1 7 



with X = El,/ Ec and C2 is small, usually negative. The function Ei (E) must 
be continuous at the boundary of the two regions. In this region, Ei has a 
much slower increasing slope on E, approaching the plateau corresponding to 
the third region, and given by the asymptotic limit. 

For particles that do not belong to the medium, but for unielement media, the 
problem can be solved in principle starting from the solution corresponding 
to particles identical with the medium ones, obtained from eqs. 5 and 6, and 
denoted by e (E). So, the equation who's solution is of interest is: 

E[ (E) ■ Se,part = J da^art [t {E - T) ~ 6 (E) + 6 (T)] (8) 



where is the electronic stopping cross section for the ion Zpart in the medium 
Zmed, and dapart is the differential cross section for an elastic scattering of the 
particle on the atom of the medium, with the stopping cross section Sn,part- In 
this situation, the division of the energy interval into three regions does not 
hold more, because, in place of Se and Sn we have now Se,part and Se,med, , 
and {Sn,part, Sn,med ) respectively. The case of media containing more than one 
element is more complicated, and is very difficult to be treated rigorously. A 
first approximation for the solution of the problem of the energy partition of a 
particle in a multielement medium is to solve separately for each component, 
and to use the average weight Bragg additivity. This method has been used 
in the case of some binary semiconductors [7,8]. 



3 Particular cases and discussion 



In the case of a medium consisting of only one atomic species, but for parti- 
cles different from the medium, the eqs. 5 and 6 could be still used, as a first 
approximation, instead of eq. 7, redefining the quantities E^ and Eh. An im- 
portant distinction is to be made between particles lighter then the medium, 
and heavier than the medium. While E^ has the same dependence on the 
characteristics of the particle and medium in both cases: 

Eh ^ 5.2 ■ 10-4 l^^-^ ^y'^J ^E^ _ (9) 



■^part K^part + ^med) 



(in MeV) E^ has different forms. 

rOr Zpart < ^nied SUd Apart < ^med- 
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E, =. 5 ■ iO'^^Ji:ill^ ^2!^ _ (10) 



Zpart [Apart + ^' 



rned) 



while for Zpart > Zmed and Apart > Amed- 



E,., - 1.25 ■ lO-^Zmed '^^"""" : ^^""^ (111 



\-^part ' -^med) 



^part 



For a given medium, consisting of only one atomic species, the family of curves 
characterising the dependence of the energy channelled into displacements, as 
a function of the recoils energy, and having the mass and charge numbers 
impart < ^med and Zpart < Z^ed) as parameters, has the following character- 
istics: 

- The maximum energy transferred into displacements corresponds to particles 
identical to the medium ones. 

- All curves start, at low energies, from the same curve; they have at low 
energies identical values of the energy spent in displacements, independent on 
the charge and mass number of the recoil, and, roughly, an dependence. 

- At higher energies, the curves start to detach from this main branch. This 
happens at lower energies if their charge and mass numbers are smaller. Then, 
the curves present a smooth increase with the energy. This means that at high 
enough energy of the incident particle, the increase of its energy determines, 
mainly, an increase of the ionisation loss. The asymptotic limit (Ep) of the 



Fig. 1. The displacement energy versus the particle (recoil) energy for C in diamond 
(continuous line), Si in silicon (dotted line), and respectively Ge in germanium 
(dashed- dotted line). The asymptotic limits of the displacement energies are also 
represented. The characteristic energies for the curves are indicated with arrows. 

equation 6 depends on the characteristics of the particle and of the medium 

as: 



Zpart impart + ^' 



med) 



In the particular case of a particle belonging to the medium, the value of 
this energy is proportional to the product of the mass and charge numbers; 
if, more, Z ~ A/2 , as is the case of C in diamond. Si in silicon and Ge in 
germanium, an A-"^ dependence is obtained for Ep. 

In Figure 1, the Lindhard curves for C in diamond. Si in silicon ad Ge in ger- 
manium are presented. The values of corresponding to the boundary between 
the first and the second energy regions for all three curves are indicated by ar- 
rows. The asymptotic limits are also shown at the highest limit of the abscissa. 
In fact, the solution given by equation 6 is the most important in the evalu- 
ation of the bulk damage produced in materials, in accelerator applications. 
For example, the standard physics programme at LHC is based on an inte- 
grated luminosity of 5-10^ pb~^, which corresponds to year of operation, for an 
annual operation time of 1.9-10^ s. The charged hadrons are produced by the 
primary interaction proton-proton at 7 GeV, while neutrons are albedo parti- 
cles. The irradiation background is continuous. The charged hadrons spectra 
(tt^, K^, p,P) simulated for various positions inside the tracking cavity of the 



Fig. 2. The displacement energy versus particle (recoil) energy for two sets of ele- 
ments near one another in the periodic table, that form two binary compounds. The 
first one indicates the particle, the second one the medium. The asymptotic limits 
(plateau) of the displacement energies are also represented. 

CMS [10] suggest that all hadrons have their kinetic energies between 10~^ - 
20 GeV. Corresponding to these kinetic energies, if only elastic interactions 
with the detector medium (diamond, silicon, germanium) are considered, the 
recoil nuclei have energies higher than some MeV, values in the saturation 
region of the Lindhard curves. The case of binary media has been treated, 
in a first approximation, as specified before, considering separately the two 
components, and then weighting with their number in the molecule. If the two 
components are adjoining each other in the periodic table, as is the case of, 
e.g., GaAs and InSb, then there is no significant difference in the Lindhard 
partition energy curves, for element placed between the two elements of the 
compound. In Figure 2, the Lindhard curves for As in Ga, Ga in As, Sb in In 
and In in Sb are given. For all of them, the asymptotic limit is also indicated. 
The curves for Ga and As could be compared with the one corresponding to 
Ge in germanium, presented in Figure 1. If the two components are far each 
other in the periodic table, then there are, in principle, two families of curves, 
one corresponding to the heaviest element, the other to the lightest one, sep- 
arated by a gap. As an illustration, in Figure 3 the results for P in In and 
In in P, as well as those for Ga in P and P in Ga are presented. These are 
useful in the evaluation of radiation defects in the semiconductor compounds 
GaP and InP. The asymptotic limit of the energy used for displacements is 
also represented. 



Fig. 3. The displacement energy versus particle (recoil) energy for two sets of ele- 
ments far one another in the periodic table, that form two binary compounds. The 
first one indicates the particle, the second one the medium. The asymptotic limits 
of the displacement energies are also represented. 

4 Summary 



Recent development of accelerator machines conduces to the necessity of the 
study of radiation effects in a multitude of materials, in intense fields of a 
large variety of particles, in a high range of energy and at different fluences. 
As essential factor in the evaluation of radiation effects is the knowledge of 
the energy partition of slowing down particles, between ionisation and dis- 
placements in the lattice of the target. The Lindhard theory gives a detailed 
description of the calculation procedure for particles identical to the medium 
and for media consisting of one atomic species, and proposes some approxi- 
mations for other situations. In the present article, we applied approximate 
methods for the calculation of the Lindhard curves for particles that do not 
belong to the medium, both for unielement and compound media. As the 
asymptotic value of these curves is of much interest for practical applications, 
the dependence of this value on the characteristics of the particle and medium 
is given. 
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